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SUfWARY 


A  study  is  made  of  the  effect  of  Initial  deviations  on  the  load 
carrying  capacity  of  thin  circular  cylindrical  shells  under  uniform 
axial  conqpression.  A  perturbation  e3q)ansion  is  used  to  reduce  the  non¬ 
linear  equations  of  von  Karman -Donnell  to  an  infinite  set  of  linear 
equations,  of  which  only  the  first  few  need  be  solved  to  obtain  a 
reasonably  accurate  solution.  The  results  for  both  infinite  shells  and 
shells  of  finite  length  indicate  that  a  small  iiiq>erf action  can  sharply 
reduce  the  maximum  load  that  a  thin -walled  cylinder  will  sustain.  In 
addition,  for  a  particular  set  of  boundary  conditions,  it  is  shown  that 
the  effect  of  the  length,  for  a  finite  shell,  is  small. 
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Governing  Equations 

The  ^.roblem  under  consideration  is  to  determine  the  effect  of 
initial  in5)erf actions  on  the  buckling  of  a  thin  circular  cylindrical  shell 
under  uniform  axial  compression.  Initially  the  infinite  shell  will  be 
treated, for  which  the  effects  of  boundary  conditions  can  be  disregarded. 
Then  the  discussion  will  be  broadened  to  include  the  solution  to  a  shell 
of  finite  length  with  classical  simple  supports. 

The  behavior  of  such  a  shell,  when  the  deformations  are  such  that 
non-linear  effects  are  important,  is  governed  by  the  well-known  von 
Karm^n-Donnell  Equations.  These  are  two  fourth -order,  non-linear,  coupled, 
partial  differential  equations.  The  first  of  the  equations  enforces 
equilibrium  in  the  radial  direction; 


dx*dy*  dx*^y* 


1  d^F*  ^ 

^  dx*^  ~ 

+ 

dy*^  dx*^ 


(1) 


where  w*  is  the  initial  deviation  from  the  perfect  cylindrical  form, 
w*  the  total  radial  displacement,  o  is  the  uniform  compressive  stress 
applied  in  the  axial  direction  and  positive  in  compression,  F*  is  the 
stress  function  that  assures  in-plane  equilibrium,  and  x*  and  y*  are, 
respectively,  the  axial  and  circumferential  coordinates.  The  shell  wal"’ 
thickness  is  denoted  by  t,  the  radius  by  R,  Young's  modulus  is  E, 
Poisson's  ratio  is  v,  and  D  is  the  bending  rigidity 
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niTRQDUCTIQN 


IDie  determination  of  the  buckling  and  postbuckllng  behavior  of  a 
thin  circular  cylindrical  shell  under  uniform  axial  compression  is  one 
of  the  most  challenging  problems  of  the  theory  of  elasticity.  The 
earliest  work,  that  of  Lorenz  [l],  Timoshenko  [2],  Southwell  [3],  and 
Fliigge  [4]  had  been  concerned  with  the  generation  of  stability  limits. 

The  lack  of  agreement  between  these  results  and  esqieriments  conducted 
during  the  1930 's  led  Donnell  [3]  and  von  Karman  and  Tsien  [6]  to 
investigate  non-linear  postbuckling  equilibrium.  They  deronstrated  the 
existence  of  equilibrium  positions  at  stresses  well  below  the  classical 
value  of  the  buckling  stress,  and  thus  were  able  to  offer  a  partial 
ejcplanation  of  the  low  loads  obtained  in  the  laboratory.  Since  that 
time  many  investigators  have  contributed  to  the  theory;  a  summary  of 
these  efforts  has  been  given  in  a  recent  survey  paper  by  Hoff  [?]•  It  is 
of  interest  for  the  present  study  to  take  special  note  of  two  aspects  of 
the  previous  work. 

The  inportance  of  initial  inaccuracies  in  the  shax>e  of  the  cylindri¬ 
cal  shell  has  been  recognized  for  a  long  time.  It  had  already  been 
studied  by  Fliigge  in  1932  and  Donnell  in  1934.  Then  in  1950  Donnell  and 
Wan  [8]  made  an  extensive  investigation  of  the  effects  of  inperfections. 
Unfortunately,  however,  their  results  may  be  viewed  with  some  suspicion 
for  in  using  a  minimization  principle,  they  minimized  with  respect  to  the 
displacement  coefficients  defining  the  initial  shape  as  well  as 
the  total  deflection.  Thus  the  system  whose  total  potential  was  minimized 
changed  its  initial  configuration  during  the  minimization.  The  same 
error  is  found  in  articles  of  Loo  [93#  Lee  [10],  and  Sobey  [ll]. 
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This  error  was  avoided  In  an  Investigation  hy  Madsen  and  Hoff  [12]. 
They  clearly  showed  that  small  Imperfections  can  have  a  significant 
effect  on  the  maximum  load  that  the  shell  can  carry.  This  maximum  load. 
Incidentally,  Is  the  only  quantity  that  can  be  observed  directly  In  a 
test  In  the  laboratory.  Further  calculations  by  Almroth  [13]  and 
Kblter  [l4],  for  axlsymmetrlc  ln5>erfectlons,  show  similar  features. 
Kolter's  i>aper  Is  a  continuation  of  an  Ingenious  thesis  [I5]  printed  In 
1945  In  which  the  principal  features  of  postbuckllng  behavior  were  derived 
from  the  properties  of  the  state  of  the  shell  at  the  crltleril  point  of 
the  small-dlsplacement  theory.  This  long-ignored  work.  In  which  the 
significance  of  the  Initial  deviations  was  fully  demonstrated  for  the 
first  time,  became  better  known  only  after  a  summary  of  It  was  published 
In  English  [I6]  In  1963. 

It  Is  also  Interesting  to  discuss  the  methods  employed  by  previous 
Investigators  to  obtain  their  solutions.  With  only  a  few  exceptions, 
beginning  with  the  early  Investigations  of  Donnell  and  von  Karman  and 
Tslen,  one  of  the  variational  principles  of  mechanics  lias  been  used, 
generally  that  of  the  Minimum  of  the  Total  Potential.  The  most  notable 
of  the  exceptions  Is  Kolter’s  analysis  of  1945  [15]*  More  recently, 

Mayers  and  Rehfleld  [17]  used  the  Relssner  variational  principle,  while 
Almroth  [13]  linearized  the  von  Karraan-Donnell  equations  and  reduced  them 
to  ordinary  differential  equations  before  Integrating  the  final  equations 
numerically  on  a  high  speed  conq)uter.  Kolter  [l4]  used  a  Galerkln 
procedure  on  the  equation  of  equilibrium  after  satisfying  the  conpati- 
bility  condition.  TTiis  is  similar  to  a  Rayleigh-Ritz  procedure.  How¬ 
ever,  only  one  term  of  the  displacement  series  was  obtained.  It  is 


vortfaHhlle  to  note  that  even  when  a  technique  as  powerful  as  the 
Minimum  of  the  Total  Potential  Ei^rgy  Is  used^  a  great  deal  of  numerical 
work  on  a  high  speed  digital  conq)Uter  Is  generally  required  to  obtain 
numerical  results. 

In  this  paper  the  effect  of  laq^rfectlons  on  the  buckling  load  of 
a  shell  will  be  studied  with  the  aid  of  a  perturbation  technique.  An 
approach  of  this  type  has  enjoysd  great  prominence  In  other  areas  of 
mathematical  physics^  especially  in  fluid  mechanics  [l8].  Few  attenq>ts 
have  been  made  to  use  this  technique  In  stability  problems.  Stein  [I9] 
has  Investigated  the  postbuckllng  behavior  of  siiiq>ly-siq>ported  elastic 
plates^  and  Kaplan  and  Fung  [20]  have  analysed  the  buckling  of  shallow 
spherical  caps^  using  perturbation  techniques.  In  both  cases  non»linear 
systems  of  equations  were  reduced  to  a  tractable  form  so  that  results 
could  be  obtained  in  a  simpler  fashion  than  previously.  In  addition^ 
good  agreement  has  been  reached  with  only  one  or  two  terms  of  the 
ejqpansion. 

In  the  present  problem  a  perturbation  expansion,  using  as  the 
perturbation  parameter  a  characteristic  amplitude  of  the  initial  deflec¬ 
tion,  has  been  used  to  reduce  the  nonlinear  von  KaiTnan-Donnell  equations 
to  an  infinite  set  of  linear  partial  differential  equations.  Only  the 
first  few  of  these  need  to  be  solved  in  order  to  obtain  a  reasonably 
accurate  solution.  Investigations  are  made  for  the  infinite  shell  and 
for  a  shell  of  finite  length  with  classical  simple  supports. 

For  the  infinite  shell,  closed  form  solutions  are  given  for  the 
deflection  function  up  to  the  third  order  of  approximation.  For  the 
finite  shell  the  satisfaction  of  the  boundary  conditions  complicates  the 
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problem  somewhat,  and  after  two  approximations  it  becomes  necessary  to 
resort  to  a  con^uter  to  obtain  terms  of  a  hijpier  order. 

It  is  shown  that  the  effect  of  Inqperfectlons  is  significant,  and 
goes  a  long  way  towards  e:q>lalning  the  lower  buckling  loads  obtained  by 
many  experimentalists.  Indications  are  also  presented  of  the  catastrophic 
failure  of  relatively  well  made  shells  which  may  buckle  at  high  loads. 
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OSEORY 


Governing  Equations 

The  problem  under  consideration  Is  to  determine  the  effect  of 
Initial  ln5)erf actions  on  the  buckling  of  a  thin  circular  cylindrical  shell 
under  uniform  axial  con^resslon.  Initially  the  Infinite  ^ell  vlU  be 
treated,  for  vhlch  the  effects  of  boundary  conditions  can  be  dlsregaz^ed. 
Bien  the  discussion  will  be  broadened  to  Include  the  solution  to  a  shell 
of  finite  length  with  classical  sln5)le  supports. 

The  behavior  of  such  a  shell,  vhen  the  deformations  are  such  that 
non-linear  effects  are  important,  is  governed  by  the  vell-knovn  von 
Karm^n-Donnell  Equations.  These  are  tvo  fourth-order,  non-linear,  coupled, 
partial  differential  equations.  The  first  of  the  equations  esaforces 
equilibrium  in  the  radial  direction: 


(D/t)  . 

a V  a V  . 

ax*®  aj’*®  ’ 


;v2  ^ 

.  O  W 


1  _ 

by*^  bi^^ 


(1) 


where  w*  is  the  initial  deviation  from  the  perfect  cylindrical  form, 
w*  the  total  radial  displacement,  a  is  the  uniform  compressive  stress 
applied  in  the  axial  direction  and  positive  in  compression,  F*  is  the 
stress  function  that  assures  in-plane  equilibrium,  and  x*  and  y*  are, 
respectively,  the  axial  and  circumferential  coordinates.  The  shell  wal"' 
thickness  is  denoted  by  t,  the  radius  by  R,  Young's  modulus  is  E, 
Poisson's  ratio  is  v,  and  D  is  the  bending  rigidity 
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The  second  equation  is  a  conqwtibility  relationship  that  ensures  the 
continuity  of  the  shell  displacements: 


vV  +  (e/r)  ^  -  (e/r)  ^  - 

bx^  Lbx*ay* 


\2  b^v*  b^v* 
bx*^by^ 


'b^v* 


^  “O  ^  *0 


ax*^  dy^ 


by*^  J 


(2) 


TOirough  the  introduction  of  suitable  dimensionless  quantities^  these 
two  equations  can  be  recast  in  the  following  form: 

4  4 

V  w  +  2p  w,  -  F,  =  ^  w»  +  q)[F,  w,  -  2F,  W;  + 

XX  'xx  0  ^  'xx  yy  xy  xy 
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„4  r  2 
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and  subscripts  after  a  comna  denote  differentiation.  It  is  easily  seen 
that  the  only  free  pairameter  reaainlng  In  the  equations  Is  p,  the 
ratio  of  the  buckling  stress  to  the  classical  critical  stress. 

The  relationship  between  the  applied  stress  and  the  average  shorten¬ 
ing  of  the  shell  In  the  axial  direction  will  next  be  given.  If  an 
average  unit  end -shortening  €  is  defined  by 


* 


then  by  means  of  the  midsurface  strain -displacement  relationship  In  the 
X*  direction 


bx*  ^  bx* 


and  the  stress -strain  law 


dS’'* 


-  -  £  + 


the  load-shortening  relationship  can  be  written, 
of  the  dimensionless  quantities  defined  earlier, 


with  the  introduction 
as 


+  |w2^^]axdy.  (1.) 


where 
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This  is  recast  once  more  into  a  slightly  more  convenient  form: 


-  X  +  f  2* 

’'A  y  Jx.  "  “ 


+  I  WQ^^]dxdy,(4a) 


where 


p  =  \/3(l-v®)'  ^  =  c/a^^ 

X  =  v/sd-v®)'  ^  =  €/(o^j^/E)  =  eA_.^  .  (4b) 


Method  of  Solution 

A  perturbation  expansion  will  be  used  to  reduce  the  nonlinear 
equations  given  above  to  an  infinite  system  of  linear  equation)^,  only  the 
first  few  of  which  need  to  be  solved  to  obtain  a  reasonably  accurate 
solution.  Since  the  primary  interest  is  in  determining  the  effect  of 
imperfections  on  the  load-carrying  capacity  of  a  shell,  it  seems  natural 
to  choose  as  a  perturbation  parameter  some  quantity  which  would  be 
indicative  of  the  degree  to  which  the  shell’s  form  deviates  from  the 
perfect  circular  cylinder.  With  this  in  mind,  the  initial  deviation  was 
taken  as 

Wq  =  e  f(x,y)  * 

where  f(x,y)  represents  the  shape  of  the  initial  deformation,  and  e 
is  a  characteristic  amplitude.  As  the  displacements  have  been  nondimen- 
sionalized  with  respect  to  the  shell  thickness,  the  number  e  actually 
represents  the  deviation  amplitude  as  a  fraction  (or  multiple)  of  the 
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thickness.  There  is  one  restriction  placed  on  the  form  function, 
f(3c,y),  i.e.,  it  must  satisfy  appropriate  boundary  conditions  in  those 
Investigations  where  shells  of  finite  length  are  examined. 

Ihen  the  characteristic  amplitude  e  is  taken  as  the  perturbation 
parameter  and  it  is  assumed  that  the  radial  displacement  w  and  the 
stress  function  F  can  be  represented  by  the  following  power  series 
In  e: 


w 


w^(x,y)e” 


(6a) 


F 


ifci 


(6b) 


It  would  seem  then  that  all  that  would  be  needed  to  effect  the 
desired  solution  would  be  to  substitute  the  expansions  (6)  into  the 
differential  equations  (3)  and  carry  forward  the  linearization  procedure. 
Then  the  first  few  terms  of  the  expansions  can  be  obtained  and  used  with 
equation  (4)  or  (4a)  to  generate  the  desired  relation  between  the  load 
and  the  end-shortening.  Unfortunately,  the  results  obtained  thusly 
would  be  spurious. 

The  curve  generated  in  this  fashion  would  have  no  true  maximum 
(curve  S,  Figure  l);  rather  it  would  approach  the  line  p  =  1  asyn^)- 
totically  as  the  shortening  parameter  X  approaches  infinity.  This 
does  not  at  all  resemble  the  type  of  behavior  expected  for  a  circular 
cylindrical  shell,  with  imperfections,  under  axial  compression. 

Fortunately  this  unexpected  result  can  be  both  explained  and 
remedied.  What  the  expansion  has  accomplished  in  the  above  approach  has 
been  the  movement  of  the  singular  point  of  the  load-shortening  curve 
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from  the  x^glon  around  p  =  X  =  1  in  such  a  way  that  now  X  becomes 
infinite  as  p  approaches  unity,  ^e  singular  nature  of  the  load- 
shortening  at  p  *  X  =  1  is  evident  from  the  cusp  in  the 
curve  (for  the  perfect  shell)  at  that  point  (Figure  l).  The  correction 
of  the  above  fallacious  results  would  seem  to  lie,  then,  in  expanding 
the  load  parameter  p  about  the  shortening  parameter  X;  and  Indeed,  a 
natural  fom  of  this  expansion  is  available  from  the  load  shortening 
relationship  (4a). 

If  the  expansions  (6)  for  the  displacement  w  and  the  stress 
function  F  are  substituted  into  the  functional  on  the  ri^t-hand  side 
of  (4a),  the  load- shortening  relation  can  be  written  in  the  following 
way; 

p  =  X  *  el,  *  .  .  .  .  .  (7) 


where 


\fl 


5(l-v^) 

A  A 
X  y 


Q  \  A  +A 
2^  ry 


riK 


yy 


V 


(7a) 


0  -A 


^2  = 


5(l-v  ) 

A  A 
X  y 


I  J^\yy  ■  ^  ''2,xx  '  I  *  I 


0  -A 


^3  ^ 


^  0  -A 


\  Av 

TT~J  J  \ 

X  y  U 

0  -A 


yy 


''4, XX  -  I  4,x  ■ 


and  so  on,  where  A  =  A  /2. 
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ll^cn  the  expansion  (?)  will  he  used  in  conjunction  with  the  expan¬ 
sions  (6)  to  linearize  the  nonlinear  differential  equations  (3). 

It  is  also  possible,  and  enlightening,  to  argue  the  use  of  the 
expansion  (7)  from  two  other  points  of  view.  The  final  aim  of  any  such 
analysis  of  the  nonlinear  behavior  of  the  circular  cylindrical  shell 
under  uniform  axial  compression  is  usually  the  determination  of  the  load¬ 
shortening  relationship  for  the  shell.  As  can  be  seen  here  from  equation 
(7)^  this  relationship  displays  the  load  parameter  p  as  a  power  series 
in  the  imperfection  amplitude  e.  Thus  it  would  seem  consistent  to  replace 
p  in  the  differential  equation  with  its  corresponding  power  series,  in 
much  the  same  fashion  as  was  done  for  the  stress  function  F  and  the 
displacement  w  by  their  expansions  (6). 

Further,  the  behavior  expected  is  illustrated  by  the  dashed  lines 
in  Figure  1.  It  is  seen  here  that  the  maxima  sought  are  relatively  close 
to  the  strai^t  line  p  =  X.  Thus  it  would  seem  only  natural  when  linear¬ 
izing  the  governing  differential  equations  to  perturb  about  this  line  at 
the  same  time. 

Both  of  these  arguments  are  helpful  for  an  understanding  of  the  prob¬ 
lem  and  its  solution;  it  should  be  remembered,  however,  that  tht  primary 
purpc'e  of  the  expansion  for  p  is  to  keep  the  appropriate  mathematical 
behavior  in  the  region  where  it  belongs. 

The  expansions  given  can  now  he  substituted  in  the  differential 
equations  to  obtain,  as  mentioned  earlier,  the  equivalent  infinite  set 
of  linearized  equations.  The  first  few  of  these  are  presented  here: 
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0(e)  \ 


V  F,  +  w,  =  f , 

1  1,XX  XX 


2,  „k 


0(e  )  O  «2  +  2X.2^^  -  =  ‘^Vl.xx  ^  ’’f^,*x’'l,3ry 


-  2F  w  +  P  w  1 
l,xy  l,xy  l,yy  IfXx-* 


o^Fg  +  ,g  «  =  -  4,, 


''o,xx“o,yy^ 


0(e^  -  Fj^^  =  -  aij^.g^^  -  2ig«,^„ 


^  ^^^l,xx^2,yy  *  ^2,xx'*l,yy  ”  ^l,xy*2,xy 


^2,xy'^l,xy  *  ^l,yy^2,xx  ^  ^2,yy”l,xx^ 

(8( 

^^3  '*’  ''3,xx  “  ^^^l,xy''2,xy  "  ''l,xx'^2,yy  "  ''2,xx''l,yy^' 


This  sequence  of  equations  may  be  written  in  the  form; 

0(e”)  +  2Xw  -  F  =  G  (x,y) 

'  n  n,xx  n,xx  n'  ' 


V  F  +  w  =  H  (x,y)  .  (9b) 

n  n^xx  n  ^  ' 

where  G^(x,y)  and  H^(x,y)  are  made  up  of  products  and  squares  of  the 
preceding  (n-l)  teiros.  It  is  in  these  expressions  that  the  nonlinearity 
of  the  original  system  is  preserved,  although  as  a  higher  order  effect 
(as  it  should  properly  be  for  this  investigation),  in  a  form  that  is 


tractable. 
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r- 


Tbe  pair  of  equations  (9)  can  be  uncoupled,  l.e.>  the  stress 
function  F^(x,y)  can  be  eliminated  between  these  two  equations  to 
yield  a  single  partial  differential  equation  in  terms  of  the  radial  dis¬ 
placement  w^(x,y) : 


+  2X 


n,xx 


+  w 

n,xxxx 


+  H 

n,xx  * 


This  equation  is  reminiscent  of  the  linearized  equation  of  Donnell 
(which  can  be  obtained  from  the  nonlinear  von  Karmw-Donnell  equations  by 
omitting  the  nonlinear  terms),  but  now  the  load  parameter  p  has  been 
replaced  with  the  shortening  parameter  X;  esad  the  equation  is  inhomo¬ 
geneous. 


The  Question  of  Covergence 

The  mathematical  foundations  of  perturbation  theory  are  not  yet 
con5)letely  established.  There  is,  at  this  time,  no  general  technique 
or  theorem  to  indicate  the  range  of  validity  of  an  asyTiq)totic  or  pertur¬ 
bation  eaqpansion.  Thus  the  engineer  must,  i.  general,  rely  on  physical 
evidence,  in  the  form  of,  say,  experimental  results,  or  different 
analytical  approaches,  in  order  to  ascertain  the  validity  of  his  solu¬ 
tions. 

In  addition,  the  following  three  indicators  seem  to  have  relevant 
implications,  both  from  the  mathematical  and  the  phj'-sical  points  of  view. 

The  first  of  these  centers  about  the  expansion  for  the  radial  dis¬ 
placement  w : 

^  2  ^  ^  n 

w  =  ew-  +  e  w»  +  .  .  .  +  e  w 


In  &  manner  similar  to  the  ratio  test  used  In  examining  the  conrargence 
of  power  series,  one  can  examine  the  behavior  of  the  ratio  (e  w  /w  .) 

In  the  neighborhood  of  the  maximum  load.  Ibe  use  of  this  Indicator, 
however,  should  be  tenqpered  by  the  knowledge  that  Immediately  after  the 
maximum  has  been  reached,  the  shell  buckles  and  the  displacements  become 
very  large.  !Hiu8  this  Indicator  Is  useful  In  the  region  just  before  the 
maximum  load. 

The  second  indicator  hinges  on  the  earlier  argument  that  the  pertur¬ 
bation  takes  place  around  the  line  p  =  X  at  the  same  time  when  the 
governing  equations  are  linearized.  Thus  an  Important  examination  would 
be  to  see  whether 

P  35  X 

near  the  maximum  load. 

The  final  test  has  a  more  physical  appeal  as  it  consists  of  an 
examination  of  the  potential  energy  of  the  system  as  given  by  each 
additional  approximation.  For  the  energy  too  may  be  written  in  the  form 

U  =  +  eU-  +  .  .  .  e“u 

0  1  n 

and  it  is  reasonable  that  if  the  perturbing  effect  is  truly  small,  the 
increase  in  energy  because  of  it  should  be  correspondingly  small. 
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SQIXmORS 


Infinite  SheU 

The  purpose  of  this  section  is  to  illustrate  the  procedure, 
described  above,  for  obtaining  the  maximum  loads  of  inperfect  cylindrical 
shells.  The  expression  "infinite  shell"  is  used  to  indicate  the  fact 
that  no  attenpt  is  made  to  satisfy  any  boundary  conditions,  and  thus  only 
particular  integrals  are  obtained  for  the  differential  equations  involved. 
Since  displacements  and  stresses  are  represented  as  trigonometric  func¬ 
tions,  vhich  automatically  satisfy  requirements  of  boundedness  at 
infinity,  this  is  a  considerable  sinplification. 

In  addition,  in  order  to  cany  out  the  integrations  indicated  in 

equations  (7)^  the  quantities  X  ,>k  are  assumed  to  represent  the  axial 

X  y 

and  circumferential  vave  lengths  of  one  buckle. 

The  initial  deviation  is  taken  in  the  form 


Vq  =  ^  ®20 

0  0 

vhere  0!^^,  are  displacement  coefficients  that  may  be  arbitrarily 

assigned,  and  m,n  are  reduced  vave  numbers. 

The  first  approximation  to  the  total  displacement  is  obtained  as 
the  solution  of 


!■  2X  ^  V-  +  V-  =  ^f  +  f , 

1,XX  1,XXXX  '3 


(11) 


A  solution  is  assumed  in  the  form 


Vi  =  cos  mx  cos  ny  +  cos  2rax  , 


(12) 
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vhich  leads  to 


(13) 

(13a) 

(13b) 


Buckling  Is  likely  to  take  place  under  the  smallest  load  if  the 
initial  deviation  pattern  is  so  selected  as  to  make  the  two  denominators 
in  equations  (l3)  vanish  simultaneously^  as  was  shown  by  Madsen  and 


Hoff  [12].  Thus  the  requirement  is 

Xii  =  X20  =  1  (I4a) 

which  leads  to 

m  =  n  =  1/2  .  (l^b) 

Then  the  initial  deviation  form  is 

f(x,y)  =  cos  |  cos  |  +  cos  x  (15) 

and  the  first  approximation  for  the  radial  displacement  is 

Wi  =  ^(“^1  cos  I  cos  I  +  a^o  cos  x)  .  (16) 

where 

i  =  (1-X)‘^  •  (16a) 
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i  . 


F-  =  -  w.  +  f, 

1  1,XX  XX 


(1?) 


After  substitution  from  (15)^  (l6)  one  finds  easily  that 


=-(l-|)f(x,y)  . 


(18) 


Now  it  is  possible  to  calculate  from  equations  (7)  the  first  approxi¬ 
mation  to  the  load-shortening  relation.  Equation  (l8)  indicates  that 
F^  is  periodic,  and  so  indeed  will  be  all  the  succeeding  approximations 
to  the  stress  function.  In  view  of  this,  the  integrals  of  (j)  can  be 
recast  as: 


ii  =  o 


(19a) 


^2- 


cp\/3(l-v^)^  ^ 


2A  A 
X  y 


f  - 


f,^)dxdy 


0  0 


(19b) 


^3  = 


r^h 

■  \\  J  J 
*  ^  0  0 


dxdy 


(19c) 


cp\A(i-v^)^  2  „  , 

■  J  J  ‘"a.x  " 

^  ^  0  0 


(I9d) 


Then  the  first  approximation  to  the  load-shortening  equation  is  obtained 


as 


oR  _  eR  j/j. 
Et  " 


eR 

t 


3(l-v^)e^ 
16 


X(2-x) 

(i-x) 


2  ^ 


(20) 


The  corresponding  stress  function  F^  is  found  from 
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One  can  proceed  in  a  similar  fashion  to  calculate  further  ajjproxi- 
mations.  In  the  present  analysis  results  were  obtained  up  to  the  third 
order  e^proximation  for  the  displacement,  and  the  results  are  given 
below: 


cos  I  cos  I  +  cos  cos  I  +  D^q  cos  x  +  cos  y^(22a) 


where 

Dll  =  -  y/sd-v^)'  -  1) 

^31  ^  25  ^31  ”  n  ^  ^^11*^20^^^  " 

^20  =  ^20  - 

^02-  -h  -  1)  . 


(22b) 

(22c) 

(22d) 

(22e) 


For  the  second  order  approximation  in  w,  the  load-shortening  relation 
becomes : 
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(23) 


^  ^  1  7,-  2x'  r  2  x(2-x)  r(„0  v2  .,  0  .2, 

Et  "  t  "  V3{l-v  )  ^ 


^  f“u  '’ll  "  =20^ 


')■ 


The  third  order  displacement  function  has  nine  harmonics: 


cos  ^  cos  I  +  G^q  cos  X  +  G^^  cos  y  +  G^^  cos  x  cos  y 


+  G^2  §■  ^  +  Gj^  cos  ~  cos  ^  +  Gj^q  cos  2x  + 


+  G,  cos  2x  cos  y 


+  G_^  cos  cos  ^  , 


(2l^a) 


where 


®11  "  ^2^  °^11  2  V3(l-v  )  ^02  “ii  ^20  ^“11^2 


+  |a3_iC2Q  +  «20^31  °^20^11  ^°^20^31  ^“20^11^ 


(24b) 


^20  ^2^  °^20  i^V3(l-v  )  I  +  “lAl  ^“11^11 


(24c) 


Go2  =  i  v/sd-v") 


(24d) 


^22  (8"' +'91!  '■“llSl  ^20^02  ?  ^°'ll^31  '*'  ^^20^2^ 


^13  (50+57^1)  ^  ■*■ 


(24f) 
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Si  “  (450+256^5  ^^^1^20  ^  ^^Sl^20  ^^2oSl  ^^“2oSl^ 


:lB?iir^  f^?iSi  KiSi^ 


(24h) 


S2  =  ilc^hu]  “iiSi  ■"  I 


(24i) 


’M^{1296)^  “^Sl  "  r  4o^51^  • 


where 


(24  j) 


I2  =  -  ^  isd-v^)  (§^-l)[(a°^)'^  + 


(24k) 


«ij  =  5°ij  - 


(s'*!) 


The  load-shortening  relationship  corresponding  to  this  order  of  the 


displacement  is 


oR  _  eR  V5(l- 
Et  "  t  TS 


P  {•’  SS 


*  (1-X)  *  ®“20°20^  ®  *  ^31  * 

"  ®4  "  2  «°  Al  "  5^  “?o«2oj)  <25 ) 

Thus  expressions  have  heen  derived  for  the  displacement  up  to  the 
order  of  0(e  ),  and  the  load -shortening  equation  to  0(e  ).  The  results 
are  presented  in  Figures  2  through  4  and  Tables  1  through  5-  The  calcu¬ 
lations  were  carried  out  on  the  Burroughs  B-5500  Digital  Computer  at  the 
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-  0  0 
Confutation  Center  of  Stanford  University.  Values  of  e,  O,^  were 

chosen  for  use  with  equations  (20),  (23)  and  (25).  Hien  the  end¬ 
shortening  parameter  X  was  stepped  from  zero  to  unity  In  intervals  of 
0.01,  and  load-shortening  curves  (Fig.  2)  were  thus  generated.  Maxima 
were  easily  obtainable  from  these  plots.  The  results  for  axisymmetric 
buckling  were  obtained  by  taking  “ll  “ 

Hie  form  of  the  load-shortening  curve  is  shown  in  Fig.  2.  It  should 
be  pointed  out  that  the  form  shown  is  not  reliable  beyond  the  maximum 
load,  because  beyond  this  point  the  validity  of  the  ejfansions  becomes 
doubtful.  Indeed,  an  examination  of  Table  3  shows  that  the  ratio 

ew  /w  - 
n'  n-1 

is  always  small  rifidit  up  to  the  vicinity  of  the  maximum,  and  then  becomes 
equal  to  or  greater  than  one  at  the  maximum  load  for  the  smaller  inferfec- 
tions.  This  latter  occurence  is  not  surprising  for  it  can  also  be  seen 
from  this  t-ble  that  the  displacement  w  increases  much  more  rapidly 
for  the  smaller  inferfections .  Ihis  implies  the  possibility  of  catas¬ 
trophic  failures,  in  the  laboratory,  of  even  well-made  specimens. 

Table  1  and  Fig.  3  show  the  variation  of  the  maximum  load  obtained 
for  a  given  value  of  the  imperfection  parameter  e  with  the  number  of 
approximations  kept  in  the  load-shortening  expression,  as  well  as  show¬ 
ing  the  variation  of  the  maximum  load  with  the  imperfection  parameter. 

From  the  viewi)oint  of  "convergence”  considerations,  there  seems  to  be 
very  good  agreement  for  small  values  of  e,  but  less  agreement  as  the 
imperfection  becomes  greater  than  one-tenth  of  the  wall-thickness  in 
anq)litude.  This  conclusion  is  also  inq)licit  in  the  results  given  in 
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Table  2,  where  it  is  seen  that 


p  s  X 

for  small  values  of  e;  but  as  the  imperfection  increases,  the  maxima 
seem  to  occur  further  away  frcmi  the  straight  line  which  reflects  the 
linear  solution.  Similar  remarks  may  be  made  with  reference  to  the 
values  of  the  potential  energies  of  the  various  approximations  which  are 
given  in  Table  4.  Incidentally,  the  potential  energy  was  obtained  by 
integration  of  the  area  under  the  load-shortening  curve. 

Plots  of  maximum  load  versus  inqperfection  amplitude  are  presented 
in  Figures  3,  4,  and  5-  The  first  of  these  is  a  graphical  con5)arison 
of  the  results  for  approximations  of  different  orders  for  nonsymmetric 
buckling.  Figure  4  shows  only  the  effect  of  ingjerfections  for  axisymmetric 
buckling.  In  the  latter  case  it  should  be  noted  that  the  linearization 
involved  in  the  perturbation  technique  is  such  that  axisymmetric  initial 
deviations  lead  to  total  deformations  that  are  also  axisymmetric. 

The  last  figure  (Figure  5)  is  a  graphical  comparison  of  results  of 
Almroth  [13]  and  Koiter  [l4]  with  results  obtained  from  the  present 
solution  for  both  axisymmetric  and  non -axisymmetric  buckling  deformations. 
The  work  of  these  two  investigators  involved  consideration  of  general 
buckling  deformations,  but  axisymmetric  initial  imperfections.  From  the 
results  of  the  present  analysis  alone  (curves  C,  D),  it  appears  that 
axisymmetric  deviations  are  the  more  critical.  Further,  for  small  values 
of  the  imperfection  amplitude  e,  Koiter's  curve  (A)  is  below  the  curve 
for  general  buckling  obtained  here.  The  c'rrve  obtained  by  Almroth,  which 
he  suggests  as  a  lower  bound  even  for  cylinders  with  non-axisymmetric 
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initial  deformations,  bounds  th3  present  general  solution  from  below, 
up  to  moderately  large  values  (0.35)  of  the  imperfection  parameter  e. 

In  addition,  some  of  the  results  obtained  by  Madsen  and  Hoff  [12] 
are  presented  here: 


p 

P 

e 

Madsen  and  Hoff 

Present  Analysis 

10"^ 

0.842 

0.894 

10“^ 

0.587 

0.610 

0.5 

no  maximum 

0.151 

For  the  first  two  cases  the  difference  is  much  less  th-in  10^.  It  might 
be  pointed  out  that  although  the  initial  imperfections  are  the  same  in 
the  two  cases,  the  total  deflected  shapes  are  different  because  higher 
order  approximations  in  the  perturbation  expansion  include  considerably 
more  harmonics  than  used  by  Madsen  and  Hoff. 

As  a  further  analysis  of  the  effect  of  axisyrametric  imperfec¬ 
tions,  the  coefficient  CK^q  was  taken  to  be  zero.  The  results  are 
shown  in  Table  5*  The  term  of  0(e  )  in  the  load-shortening  expression 
vanishes  in  this  case  so  that  we  have  fewer  values  for  comparison.  How¬ 
ever,  for  both  approximations  shown,  the  buckling  loads  are  greater  than 
or  equal  to  those  for  which  an  axisymmetric  imperfection  mode  was  present. 

Finite  Shell;  SS5  Boundary  Conditions 

With  the  perturbation  scheme  already  established,  attention  is  now 
turned  to  shells  of  finite  length.  The  solutions  obtained  thus  far  for 
shells  of  finite  length  are  eigenvalue  solutions  to  the  linearized 
Donnell  equation  (see,  for  example,  Hoff  [21]  and  Hoff  and  Soong  [22]). 
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ProB  these  solutions  critical  stresses  and  vaveforms,  hut  not  actual 
amplitude^;  can  he  obtained. 

Ihe  approach  taken  here  differ r  it  the  analysis  in  the  previous 
section  in  that  hcwiogeneous  solutions  to  the  differential  equations  are 
added  to  the  pazi;icular  integrals  to  satisfy  the  boundary  conditions. 

'Bie  boundary  conditions  considered  are  those  known  as  SS3  conditions^ 
or,  alternatively,  the  classical  sii^le  supports: 

v=w,  *0  5:v=0  (26a) 

X3C  X 

at 

X  =  t  X  .  (26b) 

Here  a  is  the  axial  normal  stress  and  v  the  reduced  circumferential 
X 

displacement  that  accai^paiy  hackling.  The  quantity  X  denotes  the 
dimensionless  shell  half-length 

X  =  (l2(l-v^)]^/‘^  (27) 

/rv 

where  L  is  the  actual  (physical)  half-length  of  the  shell. 

The  second  two  conditions  can  he  recast  in  terms  of  the  stress 
function  F.  It  follows  from  the  definition  of  the  stress  function  that 

0=0  implies  F,  =  0.  (28) 

X  yy 

Itoreover,  in  the  strain-displacement  relation 
dv*  w*  .  l/8w*\2 

one  has  w=v=0  at  x=tX.  Hence 


2k 


=  0  . 


On  the  other  hand,  from  Hooke's  law: 

%  '  I  (“y  ■  • 

But  since  It  was  stipulated  that  0=0,  and  since  c  =0  at  the 

X  Yq 

ends  of  the  shell,  then 

0=0 

y 

which  in^lies 

F,  =0  (28) 

'xx  •  ' 
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0 

where  0^  is  a  displacement  coefficient  that  may  he  arbitrarily 
assigned,  and  m,n  are  reduced  wave  numbers  given  by 


m 


n*  ^ 


(32) 


and  is  the  number  of  axial  half -wave  lengths,  and  n*  the  number 
of  circumferential  buckles. 

If  solutions  are  assumed  in  the  form 


w^  =  cos  mx  cos  ny 


cos  mx  cos  ny 


they  will  satisfy  the  differential  equations  (30)  if 


“u  " 


1  -  2X 

(  2^  2x2  2 

(m  +n  J  m 

-1 

2  /  2^  2x2 

m  (m  +n  ; 

(32a) 

(32b) 


(33a) 


2xa 


11 


H  Q  U 

4 


+  1  -  2X 


m 


2 

m 


(33b) 


With  the  forms  assumed  in  equations  (32)  the  boundary  conditions  will  be 
satisfied  if  the  physical  wave  number  in  the  axial  direction  is  odd,  for 
then 


cos  mA  =  cos  =0,  m*  =  1,3,5  ••• 


(34) 


Then  the  first  approximations  satisfy  both  the  governing  equations  and 
the  boundary  conditions. 
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If  equations  (j)  are  examined  with  a  view  towards  obtaining  the 
load-shortening  relationship,  it  is  found  that 

=  0,  (55) 

again,  for  is  still  periodic.  To  obtain  must  be  found  for 

it  is  not  yet  known  whether  Fg  will  be  periodic  for  the  finite  shell. 
If  it  is,  then  Just  w^,f  are  needed  to  calculate  Ig. 

The  details  of  further  calculations  are  omitted  here  for  the  sake 
of  brevity.  The  second  approximations  to  the  stress  function,  Fg,  and 
the  displacement,  Wg,  are  given  as 

Wg  =  (A  cos  2mx  -  cosh  q^x  cosh  q^^x)  + 

+  (B  -  cosh  Pj^x  -  (JRg^  cosh  PgX  -  cosh  p^x 

-  cosh  pgX)  cos  2ny 


F^  =  (C  cos  2mx  +  cosh  q^x  +  Bg^  +  cosh  q^x) 

+  (D  +  cosh  p^x  +  Rg^  cosh  PgX  +  R^^^  cosh  p^x 


cosh  i>2x)  cos  2ny. 


(36b) 


The  constants  A,B,C,D  are  given  as 


A  = 


8(l-8m^X+l^^) 


(37a) 


(37b) 
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(37c) 


c  = 


D  = 


8A  -  <p(a^)^n^ 
32m^ 


.m.2 
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(p  =  [I2(l-v^)]^/^ 


The  unknown  constants 


(37d) 


(37e) 


\o*  ^20'  ^21'  \l 


(38) 


are  determined  by  the  requirement  that  the  coefficients  of  each  harmonic 
in  the  circumferential  direction  must  satisfy  the  boundary  conditions 
independently.  As  the  assumption  of  hyperbolic  cosines  in  the  axial 
direction  inqplies  symmetry  in  the  axial  direction,  for  each  harmonic 
coefficient  four  boundary  condtions  (equations  (29))  are  available,  and 
four  constants  have  to  be  determined  (equation  (38)). 

The  constants  in  the  arguments  of  the  hyperbolic  functions  are 


qi  =  ~  [(1  -  X)^^^  +  id  + 
^  V2 


=  -  (I-X)^/^  -  (R  +  l6n^  -  X) 


1/2 


'4' 


■4 


=  +  (l-x)^/^  -  (R  +  l6n^  - 


=  -  (1+x)^/^  -  (R  -  l6n^  +  X)^^^ 
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(Continued) 


(39) 


Pg  =  +  (1+X)^'^^  -  (R  -  l£>n  + 

R  =  +  (1  -  32Xn^  +  256n^)^/^ 

• 

The  details  of  obtaining  these  expressions  and  the  homogeneous  solutions 
to  the  linearized  coupled  Donnell  equations  are  given  In  the  i^pendix. 

Returning  again  to  equations  (?)  It  is  apparent  that  enough 
Information  has  been  obtained  to  compute  1^.  Not  quite  so  readily 
apparent  is  the  fact  that  will  be  identically  zero;  this  result 
will  now  be  proven.  The  solutions  already  obtained  are  in  the  form 

w^  =  cos  mx  cos  ny 
cos  mx  cos  ny 
Wg  =  A(x)  +  B(x)  cos  2i?y 
Fg  =  C(x)  +  D(x)  cos  2ny^ 

Substitution  of  these  into  equations  (8c)  readily  shows  that 
will  be  of  the  form 

w^  =  E(x)  cos  ny  +  F(x)  cos  Scy 

Fj  =  G(x)  cos  ny  +  H(x)  cos  3ny. 

If  the  above  expressions  are  substituted  into  equation  ijc),  and  use  is 
made  of  orthogonality  and  periodicity,  then  it  can  be  shown  that 


Then  the  load -shortening  relationship  can  be  written  as 
p  =  X  +  +  OCe*^)  ^  (40) 
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Nov 


TT 

X  y 


n 


(p. 


>yy 


''%xx  ■  I  ''l.x  *  I 


Substitution  of  the  expressions  already  obtained  for  ^ 

yields,  after  some  manipulation. 


If  ^ 

-  V  3^1-v  ) 

-  (aKg  +  bKj^)  cosh  aX  sin  bX]  , 


where 


a  =  (1  - 

>12 


b  =  -—  (1  + 

42 


and 


K,  = 


m  cos  2mX 


K, 


^  ab[sinh^  aX  +  cos^X]  ^ 


m  cos  2raX 


^  ab[sinh^  aX  +  cos'^  bX] 


2  , ,  ^2 


with 

sinh  aX  sin  bX  -  cosh  aX  cos  bX 
Kp  =  cosh  aX  cos  bX  +  Jp  sinh  aX  sin  bX  , 


(Ul) 

(42a) 

(42b) 

(43a) 

(43b) 

(43c) 

(43d) 
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where 


=  (2X^C  -  AX  -  C)  (45e) 

Jg  =  (1  -  X^)^/^(2XC  -  A)  ,  (43f ) 


and  A,  C,  a^,  have  been  previously  defined. 

Thus  the  load-end-shoxtenlng  relationship  has  been  obtained  as  a 
function  of  four  geonetrlc  parameters:  the  Imperfection  amplitude  e, 
the  reduced  axial  wave  number  m,  the  reduced  circumferential  wave 
number  n,  and  the  dimensionless  shell  half-length  X.  In  terms  of 
physical  quantities  these  last  three  parameters  are 


n 


m 


_ m*Tr 

2[12(iI7)?^ 


L 


Numerical  calculations  have  been  carried  out  to  determine  the  effect  of 
e,m,n,X  on  the  buckling  load;  the  results  are  given  in  Tables6  through  9 
and  in  Figures  6  and  7> 

The  data  presented  in  Table  6  illustrate  the  variation  of  the 
maximum  load  with  the  anplitude  of  the  imperfection  parameter.  As  would 
be  expected  from  an  examination  of  the  results  obtained  for  the  infinite 
shell,  the  load-carrying  capacity  of  the  shell  decreases  as  the  initial 
deviation  from  straightness  increases. 
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T&bles  7,  Sf  9  are  presented  to  show  the  variation  in  maximum  load 
with  the  length  to  radius  ratio,  and  with  the  number  of  circumferential 
buckles.  The  inqperfection  parameter  was  fixed  at  e  =  0.10,  and  the 
radius -to -thickness  ratio  was  taken  as  625,  while  the  half-length-to- 
radius  ratio  was  taken  as  l/2,  1,  2  in  Tables  7,  8,  9  respectively. 

Then  the  number  of  axial  half  wave-lengths(m*)  was  varied  subje'^t  only 
to  the  restrictions  that  m*  must  be  odd  for  the  satisfaction  of  boundary 
coditions,  and  that  the  reduced  wave  number  m  could  never  be  greater 
than  unity.  This  last  restriction  arises  from  the  following  condition 
between  the  axial  wave  number  m  and  the  circumferential  wave  number  n; 

=  1 .  (kk) 

m 

This  condition  assures  buckling  at  the  classical  critical  load  for  a 
perfect  shell. 

It  is  seen  from  these  three  tables  that  the  lowest  buckling  load 
for  each  occurs  with  the  maximum  possible  number  of  axial  waves.  These 
results  are  summarized  as  follows; 


l/r 

m* 

n* 

P 

0.50 

13 

13.71 

0.672 

1.00 

27 

11.34 

0.664 

2.00 

57 

5.51 

0.652 

This  small  chart  indicates  that  as  the  length-to -radius  ratio 
is  increased,  there  is  a  small  drop  in  the  buckling  load,  and  a  sharper 
drop  in  the  number  of  circumferential  buckles. 
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vs  the 


In  Figures  6  and  7  curves  of  maximum  load  ratio  p 

lufljC 

dimensionless  half-length  X  are  plotted,  for  various  values  of  m,n. 
Again  for  any  particular  value  of  m,n,  the  length  effect  is  relatively 
small.  In  addition,  as  was  noted  above,  the  load  decreases  as  m 
increases  and  n  decreases. 
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CONCLUSIONS 


This  report  has  been  concerned  with  presenting  a  somewhat  different 
approach  to  one  of  the  most  exciting  problems  of  elasticity  theory,  the 
buckling  of  thin  circular  cylindrical  shells  under  axial  con^jression. 

The  particular  advantage  in  the  perturbation  scheme  used  here  is  that 
the  tools  of  linear  mathematics  are  now  applicable  to  the  problem.  The 
major  drawback  is  that  it  is  difficult  to  pin  down  the  exact  validity 
of  all  the  numerical  results  obtainable. 

As  far  as  the  physics  of  the  problem  is  concerned,  a  few  interesting 
trends  have  been  developed.  First  it  was  shown  that,  as  might  be 
expected,  a  shell  that  is  manufactured  with  smaller  geometric  in5)erfec- 
tions  will  be  able  to  support  more  load.  However,  when  loss  of  stability 
does  occur  for  a  nearly  perfect  shell,  the  failure  is  more  likely  to  be 
catastrophic,  for  it  has  been  shown  that  the  increase  in  the  displacements 
is  much  more  rapid  for  shells  with  small  eccentricities.  In  addition 
it  was  demonstrated  that  axisyrametric  imperfections  seem  to  be  the  most 
dangerous  to  the  load  carrying  capability  of  the  shell. 

For  shells  of  finite  length,  it  appears  that  there  is  no  appreciable 
length  effect.  The  load  does  drop  somewhat  as  the  length  of  the  shell 
increases,  but  the  only  rapid  change  occurs  in  the  number  of  circumfer¬ 
ential  buckles.  It  would  seem  that  the  longer  the  shell,  the  smaller 


the  number  of  such  buckles. 
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APPENDIX 


SOLUTIONS  TO  THE  COUPI£D,  HOMOGENEOUS  DONNELL  EQUATICaiS 


The  purpose  of  this  ^pendix  is  to  outline  briefly  the  determination 
of  the  characteristic  roots  and  the  form  of  the  solutions  to  the  coupled, 
homogeneous,  Donnell  equations: 


V  w  +  2Xw,  -  P,  =0 

XX  XX 


V  F  +  w,  =0 

XX 


Since  only  those  solutions  are  of  interest  that  are  symmetric  about 
the  centerline  of  the  shell  (x  =  O),  it  can  be  assumed  that  the  stress 
function  F  and  the  displacement  w  are  of  the  form 

F  =  R  cosh  px  cos  ny 
V  =  P  cosh  px  cos  my  ^ 

where  P,  R  are  arbitrary  constants.  Subbtitution  into  the  differential 
equations  yields  the  following  two  algebraic  equations : 

[(p^  -  +  2Xp^]P  -  p^R  =  0 

p^P  +  (p^  -  n^)^R  =  0  . 

For  the  existence  of  a  non-trivial  solution,  the  determinant  of  the  coef¬ 
ficients  of  the  above  algebraic  equations  must  vanish.  This  yields  the 
characteristic  equation  for  the  roots 

(p  -  n  )  +  2Xp  (p  -n;+p  =0 

The  eight  roots  of  this  equation  have  been  obtained  previously  in  closed 
form  in  work  done  at  Stanford  University  on  the  stability  limits  of  shells 


of  finite  length  (see,  for  exanple,  Hoff  [21]).  They  can  be  written  in 


the  form 


^  (“l  ^  Pl)'  ^2  '  ^  '“2  * 


^3  “  ^1'  ^4  “  ^3 


^5  =  '  ^6  “  ■  Pi 


P?  "  P2'  P8  "  "  Pi  • 


where 


=-(l  -  -  (R  +  4n^  - 

=  +  (1  -  -  (R  +  4n^  -  X)^/^ 

=-(l  +  X)^/^  -  (R  -  4n^  +  X)^/^ 

^2  =  +  (1  +  -  (R  -  4n^  +  X)^^^ 

R  =  +  (1  -  8n^X  +  l6n^)^/^. 

Thus  the  solutions  can  be  written  as 

w  =  ^  cosh  p^x  cos  ny 
i=l 


8 


F  ~  )  R.  cosh  p,x  c 

L  ^  ^ 

i=l 


,x  cos  ny 


As  the  hyperbolic  cosine  is  a  symnietric  function,  and  since  the  roots 
P^^Pg^Pj^Pg  are  Just  the  negatives  of  p^^Pg^Pj^Pi^^  these  solutions  can 
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actually  be  written  as 


U 

w  cosh  p^x  cci.  ny 

u 

F  =  ^  cosh  p^x  cos  ny  . 

The  boundary  conditions,  enforced  at  x  =  +  X,  are  then  automatically 
satisfied  at  x  =  -  X. 

Finally,  it  will  be  recalled  that  only  four  boundary  conditions 
should  be  stipulated,  while  it  appears  that  there  are  actually  eight 
arbitrary  constants : 


In  fact,  these  constants  are  related,  and  this  relationship  is 
determined  by  substitution  of  the  forms  of  the  solutions  for  the  displace¬ 
ment  w  and  the  stress  function  F  into  either  of  the  two  differenetial 
equations.  In  this  manner  one  obtains 


f 


where 


Q  =  -  X  +  1(1 


Thus  the  two  solutions  can  be  written  as 
w  =  -  ( QR-j^  cosh  p^x  +  cosh  pgX 


eobh  pjX  +  QR|^ 


cosn 


cos  ny 
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4 


cosh  p^x  cos  ny 


For  the  case  when  the  deformation  is  axisynnetric,  that  is  when 
n  =  0,  the  two  solutions  become 

2 

w  =  cosh  p^x  +  Pg  +  Pj  cosh  PjX  +  Pj^x 
1*  * 

The  constant  and  quadratic  terms  In  each  e3q;>resslon  are  due  to  the 
repeated  root  p  =  0.  Substitution  into  either  of  the  differential 
equations  yields: 

w  =  -  (QR^  cosh  p^x  -  ^2  ^3 

F  =  cosh  p^x  +  R^  +  Rj  cosh  p^x  . 


F  =  Rj^  cosh  p 
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TABLE 


1.  VARIATION  OF  MAXIMUM  LOAD  RATIO  P  WITH  THE  NUMBER  OF  TERMS 

max 

IN  THE  LOAD-SHORTENING  RELATIONSHIP;  INFINITE  SHELL;  a^o  = 


e 

ONE  TERM 

TWO  TERMS 

THREE  TE 

lO-** 

0.996 

0.994 

0.992 

10-5 

0.986 

0.976 

0.971 

10'^ 

0.944 

0.911 

0.894 

10-1 

0.742 

0.649 

0.610 

0.3 

0.485 

0.552 

0.314 

0.5 

0.307 

0.173 

0.151 

0.7 

0.179 

0.0637 

0.0561 

0.8 

0.129 

0.0302 

0.0284 

0.9 

0.0888 

0.00941 

0.0107 

/  0 

^0  = 

0 

mx  cos  ny  +  cos 

2iiix) 

m  = 

n  =  1/2,  = 

1.00,  =  0.25, 

u  =  0.30 

I 


TABIE  2.  COMPifflISOW  OF  RAXDfJM  LOAD  RATIO  P^ax  ^  EHD-SHOROENING 
PARAMETER  X  FOR  DIFFERENT  AFFRQXIMATIQNS  TO  LOAD¬ 
SHORTENING  RELATION. 


ONE  1 

ERM 

TWO  TERMS 
> 

THREE  TERMS 

e 

X 

P 

X 

P 

X 

P 

10-^ 

0.996-*- 

0.996 

0.994+ 

0.994 

0.993 

0.992 

10-=* 

0.990 

0.986 

0.981 

0.976 

0.973 

0.971 

10-2 

0.961 

0.944 

0.928 

0.911 

0.907 

0.894 

10-1 

0.825 

0.742 

0.745 

0.649 

0.676 

0.610 

0.3 

0.644 

0.485 

0.488 

0.352 

0.429 

0.314 

0.5 

0.495 

0.307 

0.310 

0.173 

0.264 

0.151 

0.7 

0.363 

0.179 

0.180 

0.0637 

0.149 

0.0561 

0.8 

0.314 

0.129 

0.120 

0.0302 

0.0990 

0.0284 

0.9 

0.248 

0.0888 

0.0599 

0.00941 

0.0660 

0.0107 

/  0  Ox 

Vq  =  cos  mx  cos  ny  +  cos  mx; 

m  =  n  =  1/2,  =  1.00,  0°^  =  0.25;  v  =  0.30 
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TABI£  3.  COMPARIoOM  OF  DISPLACEMENTS  AT  MAXIMUM  LOAD 


e 

X 

^2 

”1 

!!3 

“2 

V 

lo-** 

O.ST^* 

0.097 

0.193 

0.00554 

0.990 

0.828 

1.653 

0.04628 

(M)  0.993 

1.631 

3.256 

0.161 

10“^ 

O.9IH 

0.179 

0.353 

0.0267 

0.957 

0.352 

0.696 

0.0479 

(M)  0.973 

0.971 

1.929 

0.191 

10-2 

0.875 

0.382 

0.738 

0.167 

0.892 

0.512 

0.994 

0.234 

(M)  0.907 

0.720 

1.404 

0.374 

10"^ 

0.644 

0.392 

0.702 

0.586 

0.660 

0.438 

0.790 

0.658 

(M)  0.676 

0.492 

0.892 

0.748 

0.3 

0.396 

0.284 

0.438 

0.875 

0.413 

0.310 

0.486 

0.933 

(m)  0.429 

0.339 

0.538 

1.000 

0.5 

0.231 

0.183 

0.213 

0.993 

0.248 

0.203 

0.249 

1.042 

(M)  0.264 

0.225 

0.287 

1.095 

0.7 

0.116 

0.101 

0.0.334 

1.093 

0.132 

0.119 

0.0646 

1.136 

(M)  0.149 

0.139 

0.0972 

1.184 

0.8 

0.0660 

0.0605 

-0.0534 

1.132 

0.0826 

0.0778 

-0.0246 

1.173 

(M)  0.0990 

0.0962 

0.0061 

1.218 

0.9 

0.0330 

0.0321 

-0.1191 

1.196 

0.0495 

0.0496 

-0.0905 

1.237 

(M)  0.0660 

0.0680 

-0.0600 

1.282 

(M)  denotes  point  where  maximum  load  occurs. 


u  ^  at  nr 
■0  -''11 


^sjo  *jjr  ^  ^^20  e^UlXJ 


CL,  =  1.00,  a°  =  0.25,  m  =  n  =  l/2,  u  =  0.30 


11 
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TABI£  4.  COMPARISOW  OF  POEEMTIAL  EMERGY  OF  VARIOUS  APPROXI¬ 
MATIONS  UP  TO  MAXIMUM  LOADS. 


e 

X 

"(0) 

^(1) 

^(2) 

"(3) 

0 

1.000 

0.184 

0.184 

0.184 

0.184 

10-“ 

0.993 

0.181 

0.181 

0.181 

0.181 

10-2 

0.973 

0.174 

0.174 

0.174 

0.174 

10-2 

0.907 

0.151 

0.151 

0.151 

0.151 

10*^ 

0.676 

0.0840 

0.0827 

0.0820 

0.0815 

0.3 

0.429 

0.0338 

0.0311 

0.0294 

0.0288 

0.5 

0.264 

0.0128 

0.0106 

0.00910 

0.00882 

0.7 

0.149 

0.00405 

0.00286 

0.00198 

0.00195 

0.8 

0.0990 

0.00180 

0.00114 

0.000639 

0.000641 

0.9 

0.0660 

0.00080 

0.000445 

0.000152 

0.000179 

Vq  =  cos  mx  cos  ny  +  cos  2mx) 

m  =  n  =  1/2,  =  1.00,  =  0.25,  v  =  0.30 

^(m)  ^0  ®^1 


TABIiE  5.  VARIATION  OF  MAXIMUM  LOAD  RATIO  P  WITH  THE  NUMBER  OF 

max 

TERMS  IN  THE  LOAD-SHORTENING  RELATIONSHIP;  INFINITE  SHELL; 


e 

ONE  TERM 

THREE  TERMS 

10-“ 

0.996 

0.99!^ 

10-" 

0.990 

0.976 

10-2 

0.9^18 

0.914 

10-1 

0.77i* 

0.680 

0.3 

0.5i^5 

0.432 

0.5 

0.382 

0.281 

0.7 

0.258 

0.180 

0.8 

0.206 

0.143 

0.9 

0.161 

0.113 

/  0  Ox 

Wq  =  cos  mx  cos  ny  +  2mx; 

m  =  n  =  1/2,  =  1.00,  =  0.00,  u  =  0.30 


S 


! 


TABLE  6.  INFLUENCE  OF  IMPERFECTION  SIZE  ON  MAXIMUM  LOAD  RATIO  OF 
FINITE  LENSTH  SHELL 


X 

e 

P 

1.7^5 

0.10 

0.694 

0.010 

0.99+ 

0.0010 

0.99+ 

26.180 

0.10 

0.672 

0.010 

0.928 

0.0010 

0.99+ 

99.484 

0.10 

0.671 

0.010 

0.928 

0.0010 

0.99+ 

''o 

=  e  cos  mx  cos  ny 

m  =  0.900, 

0 

n  =  0.300,  =  1, 

V  =  0.30 
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TABI£  7.  VARIATION  OF  MAXIMUM  LOAD  RATIO  P  WITH  NUMBER  OF 

max 

BUCKLES;  L/R  =0.50 


m* 

n* 

m 

n 

P 

1 

11.55 

0.069 

0.254 

C.939 

3 

18.42 

0.207 

0.405 

0.872 

5 

21.61 

0.346 

0.476 

0.822 

7 

22.71 

0.484 

0.500 

0.786 

9 

22.03 

0.622 

0.485 

0.747 

11 

19.40 

0.760 

0.427 

0.708 

13 

13.71 

0.899 

0.302 

0.672 

e  =  0.10,  R/t  =  625,  L/R  =  O.5O,  X  =  22.72 


^7 


WITH  NUMBEIR  OF 


TABIE  8.  VARIiO’ICMI  OF  MAXIMUM  LOAD  RATIO  P 
BUCKLES;  L/R  =  1.00 


m* 

n* 

m 

n 

P 

1 

8.30 

0.0346 

0.183 

0.961 

3 

13.86 

0.104 

0.305 

0.920 

5 

17.18 

0.173 

0.378 

0.888 

7 

19.46 

0.242 

0.428 

0.860 

9 

21.04 

0.311 

0.463 

0.835 

11 

22.06 

0.380 

0.485 

0.812 

13 

22.61 

0.449 

0.497 

0.793 

15 

22.71 

0.518 

0.500 

0.77  V 

17 

22.37 

0.588 

0.492 

0.75V 

19 

21.58 

0.657 

0.475 

0.734 

21 

20.27 

0.726 

0.446 

0.716 

23 

18.35 

0.795 

0.404 

0.697 

25 

15.57 

0.864 

0.343 

0.680 

27 

11.34 

0.933 

0.250 

0.664 

e  =  0.10,  R/t  =  625,  L/R  =  1.00,  X  =  45.45 
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Figure  1.  Generalized  Load  vs  End  Shortening 


NON  -  SYMMETRIC  BUCKLINO 
W  -  0.30 


(A)  ONE  TERM  * 

(B)  TWO  TERM 

(C)  THREE  TERM 


Figure  3.  Maximum  Load  Ratio  vs  Imperfection  An5)litude;  Non- 

Syinmetric  Buckling 
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(B)  ALMROTH'S  LOWER  BOUND 

(C)  PRESENT  SOLUTION; 

NON-AXiSYMMETRiC 

(D)  PRESENT  SOLUTION; 

AXISYMMETRIC 


Figure  5.  Maximum  Load  Ratio  p  vs  Imperfection  Amplitude;  Comparison 

Itl&X 

with  the  Results  of  Koiter  and  Alraroth 


arsor 


Figure  7-  Maximum  Load  Ratio  p  vs  Half-Length 
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